A Reappraisal of the Electric Dipole Moment Enhancement Factor for Thallium 
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The electric dipole moment (EDM) enhancement factor of atomic Tl is of considerable interest as 
it has been used in determining the most accurate intrinsic electron EDM limit to date. However, 
the value of this quantity varies from —179 to —1041 in different approximations. In view of the 
large uncertainties associated with many of these calculations, we have employed the relativistic 
coupled-cluster theory with single and double excitations and a subset of leading triple excitations 
and obtained the EDM enhancement factor of TI as —466, which in combination with the most 
accurate measured value of Tl EDM yields 2.0 x lO^^^e cm as the new upper limit for the EDM 
of the electron. The importance of all-order correlation effects is emphasized and their trends are 
compared with those of two other ab initio calculations. 
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The EDM of a nondegeneratc physical system [l|, |2| 
arises from the sim.ultaneous violations of parity {V) 
and time-reversal (7") symmetries. The invariance of 7~- 
symmetry is also associated with the invariance of the 
CV (combined charge conjugation (C) and V symmetries) 
symmetry on the basis of the CVT theorem j3|. Thus, 
EDMs of atoms can shed hght on CV violations origi- 
nating in the leptonic, semi-leptonic and hadronic sec- 
tors [J, l5|i some of which are not well understood so far. 
The knowledge of EDMs provide valuable insights into 
some profound questions such as the existence of new 
physics beyond the standard model (SM) and the matter- 
antimatter asymmetry in the universe [^-|6[. Given the 
current interest in understanding different types of CV 
violations originating from elementary particles both by 
accelerator and non-accelerator based approaches, EDM 
searches are of great significance. 

The EDMs of paramagnetic atoms are sensitive to the 
the EDM of the electron [7| . The most accurate limit for 
the latter has been obtained by combining the results of 
the Tl EDM measurement and the EDM enhancement 
factor (EF), defined as the ratio of the EDM of the atom 
to that of the electron, of this atom. Liu and Kelly have 
calculated this EF, and found its value to be —585, with 
an error bar of 5 — 10%, using a linearized version of the 
relativistic coupled-cluster (RCC) theory [8[. A more re- 
cent calculation of this quantity for Tl by Dzuba and 
Flambaum uses a hybrid approach combining the con- 
figuration interaction (CI) method and many-body per- 
turbation theory (MBPT)and reports —582 with an es- 
timated accuracy of 3% Q. Both these results, coinci- 
dentally, are in good agreement with each other. Nev- 
ertheless, considering the fact that the values of the Tl 
EDM EF reported in the literature range from —179 to 
— 1041 [ISHITI, there is clearly a need for high precision 
calculations of this quantity. Thus, the primary focus of 



this paper is to determine the EF of Tl by proceeding 
beyond [8|, |9[. An accurate treatment of the unusually 
strong electron correlation effects in the ground state Tl 
EF warrants the use of an all-order relativistic many- 
body method like the RCC theory. 

The open-shell RCC theory with single, double and a 
subset of leading triple excitations employed in the cal- 
culation of the EDM EFs is discussed in detail in |10l - [ia | . 
However, we briefly present below the salient features of 
this method for the sake of completeness. The effective 
one electron form of the interaction Hamiltonian due to 
the electron EDM is given by [l^] , 
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where /3 and 7^ are the usual Dirac matrices, pj is the 
momentum vector of the j*'' electron. As the strength 
of the EDM interaction is sufficiently weak, we consider 
the wavefunction expansion only up to first-order in per- 
turbation. Thus, the modified atomic wavefunction for a 
valence electron '1;' is given by. 
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In the RCC theory, the unperturbed and perturbed 
wavefunctions can be expressed as 10l4l2l|. 



and 
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where |$„) is the Dirac- Fock (DF) wavefunction obtained 
by appending the valence electron v to the closed-shell 
([5(i^°] 6s^) reference wavefunction, T^^' and Sv are the 
excitation operators for core and valence electrons in an 



unperturbed case, where as, T^^^ and Si, are their first 
order corrections. The atomic wavef unctions are calcu- 
lated using the Dirac- Coulomb Hamiltonian given by, 
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where a and /3 arc Dirac matrices, Vn{ri) is the nuclear 
potential and Vc{rij) is the Coulomb potential on the 
electron i due to the j*^ electron. 

We consider only the single and double excitation oper- 
ators in the expansion of the RCC wavefunctions (termed 
as CCSD approximation), by defining. 



T^Ti+T2 and S^ = Siv + 821,, 
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for both the perturbed and unperturbed operators. Fur- 
ther, we construct triple excitation operators for Sv ' as. 
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which are used to evaluate CCSD amplitudes iteratively. 
This is referred to as CCSD(T) method. Here, e^ is the 
single particle energy of an orbital i. 

The final expression for the EDM EF (7^ = -2^) in 
terms of the coupled-cluster operators is given by, 

1 ($.„|{1 + si"^'}D(^{T(^\l + 5^°') + ^l''}|$.) 
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e^''"'De^'"\ndD = 



where the dressed operator £)(") 
ef is the electric dipole moment operator due to the 
applied electric field. The procedure for the calculation 
of the above expression is discussed elsewhere 10l4l2| . 




FIG. 1: The leading correlation diagrams: (i) DT^ , 

(ii) D SlJ , (iii) D S2J ■ The exchange and hermitian 

conjugate diagrams are not shown. Labels v, a and p 

refer to valence, core and virtual orbitals, respectively. 

In Table HI we present the contributions from different 
RCC terms along with the lowest order (DF) contribu- 
tions to TZ and compare them with the results reported 



TABLE I: The contributions from various RCC terms 

to the EDM EF (7e) of Tl. A new quantity, X ^ T + Sy 

is defined to compare our results with those of [8J. 
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Term 


This Work Liu 


& Kelly [8] 


\J-^ J- 1 jlowcst-ordcr 


-153.6 


-153.2 


\J-^ ^\v jlowcst-ordcr 


-268.5 


-267.3 


\^ J- 1 j higher-order 


-224.7 


-342.1 


\^ >^lv jhighcr-ordcr 


-45.5 


-102.5 


7?X(1' 


248.0 


240.9 


xf^'oxi^^ 


22.5 


22.4 


xf^'oxi^^ 


-78.2 


49.3 


xf^'ox^^^ 


21.5 


-36.9 


xf^'oxi^^xf^ 


-4.3 


-2.2 


Higher order RCC terms 


13.0 


— 


Normalization contribution 


3.8 


6.5 


Total EDM EF 


-466 


-585 



by Liu and Kelly [8[. It is evident from this table that 
the bulk of the contributions to TZ comes from the RCC 
terms DT[''\ D S^^J and D 8^2^ ■ These important all- 
order correlation effects involving the core, valence and 
core- valence sectors are shown diagrammatically in FiglT] 
The single largest contribution (— 378), to TZ comes from 
the all-order core correlation effect of which its lowest or- 
der contribution is less than half of its magnitude. The 
all-order valence correlation contribution from D SlJ is 
— 314 while its DF contribution is — 269. The differ- 
ence between the all-order and the lowest-order contri- 
butions of the terms D S{J and D T^ ' clearly indicate 
that the all-order core correlations are of crucial impor- 
tance in the case of Tl. It is indeed significant that one of 
the classes of the all-order core-polarization effects repre- 
sented by D S2J makes the largest positive contribution, 
of magnitude 248. This dramatically reduces the final 
result. There is a non-negligible contribution from the 
other higher order terms; however, many of them cancel 
each other to give an effective value of only about 13. 
The normalization of the RCC wavefunction also gives a 
contribution of 4 and the total result for the ground state 
EDM enhancement factor for Tl due to the intrinsic EDM 
of the electron amounts to —466. 

The H'^fij^ being an odd-parity operator mixes the 
atomic states of opposite parities, however, with the same 
angular momentum. Therefore, we have investigated the 
role of various intermediate states of s symmetry. The 
RCC terms such as: DT[^\ D S^^] , S^^^^ D S^^J and 
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S2y D S{J , have got significant contributions from the 
orbitals, from 5s through lis and are shown in Fig. [2l 
For comparison, we have also shown the DF contributions 
from these individual orbitals. The correlation effects for 
the 6s orbital are much larger than its DF contribution 



itself. The core electrons of d symmetry also influence 
the RCC amplitudes indirectly. High accuracy calcula- 
tions should therefore employ an all-order method such 
as the RCC theory to account for the strong correlation 
effects in a comprehensive manner. 
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FIG. 2: The contributions from singly excited 

intermediate states, 5s through lis, in DF and RCC 

approximations are compared. 

We have also investigated the role of different doubly 
excited states that contribute through D S2J . These con- 
tributions are given in Table |lTl As seen from the table, 
the largest contribution comes when the 6 s orbital is ex- 
cited to the valence orbital, 6pif2 followed by the latter's 
fine structure partner, the 6^3/2 orbital. There arc also, 
however, some non-negligible contributions coming from 
the 5d orbitals when they get excited to the virtual P3/2 
and /5/2,7/2 orbitals. 



TABLE II: Contributions from selected doubly excited 

(1) 



states through D 82^ term 



core(a) virtual(p) 


Results 


core(o) virtual(p) 


Results 


6s 


6pi/2 


170.88 


5d3/2 


8f5/2 


0.92 


6s 


7Pl/2 


3.61 


5d3/2 


9f5/2 


2.06 


6s 


8Pl/2 


2.01 


5d5/2 


6P3/2 


5.37 


6s 


9pi/2 


4.50 


5d5/2 


7p3/2 


0.76 


6s 


6P3/2 


31.42 


5d5/2 


8P3/2 


0.59 


6s 


7p3/2 


3.00 


5d5/2 


9p3/2 


1.46 


6s 


8P3/2 


1.95 


5d5/2 


8f7/2 


2.34 


6s 


9p3/2 


5.26 


5d5/2 


9f7/2 


5.36 



The EDM enhancement factor result for Tl from dif- 
ferent calculations are compared in Table ITlIl It is quite 
apparent that the published results lie in a rather large 
range. Although, the overall trends of the majority of the 
correlation terms in [8| appear similar to ours, there is a 
significant difference in the magnitude of many of those 



correlation contributions. This can be due to the several 
approximations considered in their calculation, some of 
which are the following: (i) they have considered only 
the one-body form for the unperturbed atomic Hamilto- 
nian and as well as for the EDM Hamiltonian thereby ne- 
glecting the important contributions partly from the DF 
potential and largely from the two electron Coulomb in- 
teraction, (ii) they have considered only the linear terms 
and a few selected non-linear terms in their calculations, 
(iii) they have included selective triple excitations ap- 
proximately in the unperturbed singles amplitude equa- 



.(0) 



(0) 



tions only (i.e., in T^ and S\J equations), where as, 
they are completely ignored in the unperturbed doubles 
equations, (iv) they have frozen the inner core up to the 
4s orbital for the calculation of unperturbed amplitudes, 
where as, for solving the perturbed doubles equations 
they further freeze 4s, 4p and 4d orbitals. In contrast to 
the above set backs, we consider all the non-linear terms 
arising from the single and double excitations. We also 
consider the leading triple excitations in both the unper- 
turbed singles and doubles cluster equations. We solve 
both the unperturbed and perturbed, closed- and open- 
shell equations, self-consistently. We have performed the 
all-clcctron relativistic CCSD(T) calculation and the to- 
tal active orbitals considered here, are: 14s, 13pi/2,3/2, 
13^3/2,5/31 9/5/2,7/2 and 837/2,9/2- 

The comparison between the results of Dzuba and 
Flambaum based on a combined CI-I-MBPT approach [9| 
and our all-order CCSD(T) is not straight forward. They 
have considered the three outer shell electrons (6s^ 6^1/2) 
as valence and the rest as core. The core, virtual and 
valence orbitals are generated in a V^-'^"^ closed-shell po- 
tential in contrast to the V'^^^ potential used in our cal- 
culation. The valence- valence correlations are evaluated 
by CI, while the valence-core and the core-core correla- 
tions by MBPT. In view of the strong correlation effects 
in Tl, it is important to treat them by an all-order corre- 
lation method instead of dealing with it as a three valence 
system within the framework of an hybrid CI and finite 
order MBPT approach. It appears from the previous 
work of Dzuba and Flambaum that the V &iT violating 
Hamiltonian used in [9J , considers only the internal elec- 
tric field due to the nucleus and not the electrons. The 
major drawback of the latter work is that all the corre- 
lations obtained by MBPT (except for s-electrons in the 
one-body correlation operator Si [9|) are considered only 
up to second order although they have emphasized the 
importance of valence-core correlations in their paper. In 
contrast, we have considered all these effects to all-orders 
in the residual Coulomb interaction in the framework of 
the fuU fledged CCSD(T) theory. The CI-hMBPT result 
of [9(] compares well with that of [8|. The latter as men- 
tioned earlier is based on a linearized CCSD(T) approach 
with several approximations. However, the agreement 
between these two results is fortuitous. 

We have also compared our result with a few ab initio 



TABLE III: Comparison of the ground state EDM EF 
(TV) of Tl from different calculations. 



7^ 



Method 



Reference 



-466(10) 

-582(20) 

-585(30 - 60) 

-179 

-301 

-502,-562,-607,-1041 MBPT(l) 

-500 semi-emp. 

-716 semi-emp. 



CCSD(T) This Work 
CI-fMBPT Dzuba 2009 [9J 
LCCSD(T) Liu 1992 [8] 
MBPT(2) Hartley 1990[13l 
MBPT(2) Kraftmakher 1988[14| 
Johnson 1986 fl5] 
Flambaum 1976 fl6] 
Sandars 1975 fl7] 



finite-order MBPT results |13l - lla | and semi-empirical es- 
timations [leJlSin Table Uni Apparently, the MBPT 
calculations 13|, \l^ have under estimated the results, 
where as, the first-order MBPT calculation with semi- 
empirical potentials [l5| and the semi-empirical calcula- 
tions 



16|, ll7| have over estimated the EDM EF of Tl. 



The results of a few relevant physical quantities such as 
the electric dipole (El) matrix element for 7S — >■ 6P1/2 
transition and the magnetic dipole hyperfine structure 
constants of the ground state 6P1/2 and the lowest singly 
excited state 75 have been presented in Table |TVl It is 
clear that, our results are in better agreement with the 
experiments than those reported in [9|. 



TABLE IV: Comparison of the El matrix element and 
hyperfine structure constants of low-lying states in Tl. 

This Work Dzuba fc Flambaum [9] Exp"t 
Transition 
7S^6Pi/2' 
State 

6P1/2 21053 ^1067 21311 I19J 

7S 11894 11417 12297 [20] 



El tr. amp. 
T82 L73 1.81(2) |18] 

Magnetic dipole hfs constant 



The error estimates quoted in Q for Tl EDM EF does 
not seem to be convincing for the reason that they have 
treated the valence-core correlations for the most part 
only up to second order in perturbation theory. Clearly 
an all-order correlation calculation is needed for a reliable 
error estimation. Although, in the CCSD(T) approach 
the effect of the dominant triple excitations are consid- 
ered while evaluating the CCSD amplitudes, their direct 
contributions arc ignored while computing Eq. ([8]). To 
estimate the maximum possible error on the EF result, 
we evaluate the contribution coming from the triple ex- 
P9-'(o) o.^ QPi-'W where the latter is 



citations, such as, 5*^^^ and S: 



defined as (for details, see [21 
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Ho 5*21, 



Considering the error due to the basis set incompleteness 
and the error due to the neglected triples, we obtain the 
final result for the EDM EF of Tl to be -466(10). By 
combining the EDM experimental result for Tl, —(4.0 ± 
4.3) X 10^^^ e cm, given by Regan et al. 12^ and our ac- 
curate theoretical EDM EF, we obtain a new limit for 
the electron EDM to be (8.7 ± 9.3) x lO-^^ecm which 
translates into an upper limit, de < 2.0 x IO^^^cctti at 
90% confidence level. In conclusion, we have obtained 
the most accurate limit for the electron EDM to date 
by improving the value of the EDM EF for the ground 
state of atomic Tl by treating the unusually strong elec- 
tron correlation effects to all-order using the relativistic 
CCSD(T) method. 
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